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We study some bulk thermodynamical characteristics, meson properties and the nucleon as a 
baryon number one soliton in a hot quark matter in the NJL model as well as in a hot nucleon 
matter in a hybrid NJL model in which the Dirac sea of quarks is combined with a Fermi sea of 
nucleons. In both cases, working in mean-field approximation, we find a chiral phase transition from 
Goldstone to Wigner phase. At finite density the chiral order parameter and the constituent quark 
mass have a non-monotonic temperature dependence - at finite temperatures not close to the critical 
one they are less affected than in the cold matter. Whereas the quark matter is rather soft against 
thermal fluctuations and the corresponding chiral phase transition is smooth, the nucleon matter 
is much stiffer and the chiral phase transition is very sharp. The thermodynamical variables show 
large discontinuities which is an indication for a first order phase transition. We solve the B = 1 
solitonic sector of the NJL model in the presence of an external hot quark and nucleon medium. In 
the hot medium at intermediate temperature the soliton is more bound and less swelled than in the 
case of a cold matter. At some critical temperature, which for the nucleon matter coincides with the 
critical temperature for the chiral phase transition, we find no more a localized solution. According 
to this model scenario one should expect a sharp phase transition from the nucleon to the quark 
matter. 

PACS numbers: 12.39.Fe, 14.20.Dh, 14.40.Aq, ll.10.Wx 

I. INTRODUCTION 

At some finite density and/or temperature it is generally expected a restoration of the chiral symmetry and a 
deconfinement, and hence a change of the structure of the hadrons immersed in a hot and dense medium. It is a topic 
of an increasing interest related to the evolution of the early Universe as well as to the processes in the interior of 
the neutron stars. Rather encouragingly, direct experimental studies of such phenomena in the relativistic heavy-ion 
reactions are also possible. Indeed, the lattice QCD calculations suggest phase transitions concerning both chiral 
symmetry restoration and color deconfinement at relatively low temperatures (of about 150 - 160 MeV) jjJI). These 
low temperatures are experimentally accessible even in the existing accelerators and the corresponding heavy-ion 
experiments are already planned at BNL and CERN || . From the side of theory the low critical values mean that we 
have to deal with non-perturbative phenomena. Since the analytical as well as the numerical (lattice QCD) methods 
are not developed enough to allow fully the solution of low-energy non-perturbative, especially if the baryons are 
involved, it is a motivation to apply effective quark models. The Nambu - Jona-Lasinio (NJL) model [Q is the 
simplest purely quark model which incorporates the chiral symmetry and allows for its spontaneous breaking but it 
lacks confinement. It shows, however, a significant success in description of meson structure as well as of the properties 
of the nucleon as a non-topological soliton. Quite encouragingly, assuming that the nuclear medium can approximately 
be replaced by a uniform quark medium, the model describes (for review see ^) a restoration of the chiral symmetry 
at finite temperature and/or density in a quantitative agreement with the lattice QCD calculations jl],^] as well as with 
the chiral perturbation theory Following a scenario in which the restoration of the chiral symmetry in medium 
is considered as a relevant mechanism for a scale change which modify the hadron structure, the NJL model seems 
to provide an appropriate working scheme to study the hadron properties in hot medium. Assuming that the baryon 
medium can be approximated by a quark medium the bulk thermodynamical variables and medium modified meson 
properties has been repeatedly studied |^|,0J|J|,|l^,|ll|,^2|,^3| . The nucleon as a B — 1 soliton in a cold quark medium 
has been already investigated pi . The results show that the nucleon as a soliton at finite density is less bound and 
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increases in size, and at some critical value of about two times nuclear matter density the nucleon as a soliton does 
not exist. 

In the present paper our task is to study some bulk thermodynamical characteristics and meson properties at finite 
temperature and density as well as the medium modification of the nucleon as a bound state of N c quarks coupled 
to the polarized medium (Fermi sea) and Dirac sea. In order to respect to some extent confinement we consider 
two different cases, namely a medium of nucleons, which is relevant before the phase transition, as well as a quark 
medium. In the former case, similar to the approach of Walecka p5| , the nucleons are coupled directly to the meson 
fields. However, the NJL model differs from the Walecka approach in the meson fields, which appear basically as 
quark-antiquarks excitations, and also in the way in which the chiral symmetry is dynamically broken. Such a model 
picture has been considered by Jaminon et al. in the case of cold matter. Solving the B — 1 solitonic sector 
of the NJL model in the presence of an external hot medium in a self-consistent way with the polarization of both 
Fermi and Dirac sea taken into account, we study the modification of the nucleon structure due to the medium. 
We work in mean-field approximation, which means that the meson quantum (loop) effects are not included in our 
considerations. Since the meson loop effects are dominant at low temperatures and vanishing density (pions are the 
lightest mode) |n]|^] , we restrict ourselves to consider mainly the case of finite density and large enough temperatures 
where the nucleon and quark degrees of freedom are most relevant. This case is also of interest for the relativistic 
heavy-ion reactions. 



II. NJL MODEL AT FINITE TEMPERATURE AND DENSITY 

We use the Nambu-Jona-Lasinio (NJL) model [Q which is the simplest purely quark model describing spontaneous 
chiral symmetry breakdown. The SU(2)-version of the NJL lagrangian contains chirally invariant local scalar and 
pseudoscalar four-quark interaction: 



G 

C = ^{ip-mo)^ + — [(*W) 2 + (^f 75 *) 2 ] 



(1) 



where \1/ is the quark field, G is the coupling constant, f are the Pauli matrices in the isospin space and toq is the 
current mass taken equal for both up and down quarks. Applying the well known bosonization procedure [[15| the 
NJL model is expressed in terms of auxiliary meson fields <t, 7?: 



£ = * ( ip-a-in -T75 ) * - ^(ct 2 + if 2 ) + . 



(2) 



Using functional integral technique the quantized theory at finite temperature and density can be written in terms 
of the corresponding euclidean grand canonical partition function [ftq] : 



Z = Tr expj-/3(# - fiN)} = J V^V^ expj^ dr J d 3 x(C - #V*)} , 



(3) 



where V is the volume of the system, (3 is the inverse temperature and \i is the chemical potential. The integration 
over the quarks can be done exactly, whereas for the integration over the mesons we use a large N c saddle-point 
(mean-field) approximation. It means that the meson fields are treated classically - no meson loops are taken into 
account. Following refs. [l6| ] we replace the integration over the imaginary time by a sum over fermionic Matsubara 
frequencies k — > (2n + l)ir/[3. Finally we get for the effective action 
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The energies e a are eigenvalues of the one-particle hamiltonian h: 
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and $„ are eigenfunctions. The saddle-point solution makes the effective action stationary 
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with the number of particles N in the volume V 



N = _ldS eff 
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kept fixed. Here a c and tt c are the "classical" values of the meson fields and /i is the chemical potential related to the 
number of particles N . 

The thermodynamical characteristics of a many-body system are specified by the thermodynamical potential 

_ S e ff(ji,P) 
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It should be noticed that the saddle-point solution (a Cl n c ) minimizes not Q but the Hclmholtz free energy 



F = ft - n 



(8) 
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with a constraint (0) and /i playing a role of a Lagrange multiplier. 

In the mean-field approximation (leading order in 1/N e ) the inverse meson propagator is given by the second 
variation of the effective action at the stationary point cr c , 7? c : 
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where (j) stands for both meson fields. The on-shell meson masses correspond to the poles of the meson propagator 
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at q = and the physical quark meson coupling constants are given by the residue of the propagator at the pole 
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lim y+mDK^q 2 ). 



(12) 



Due to the local four-fermion interaction the lagrangian (Q) is not renormalizable and a regularization procedure 
with an appropriate cut-off A is needed to make the effective action finite. Actually only the part of the effective 
action S e ff(/i = 0, = oo), coming from the Dirac sea (negative-energy part of the spectrum), is divergent and we 
regularize it using the proper-time regularization 
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It is easy to check that the difference S e ff(/i,P) — S e ff(fi — 0,f3 = oo) is finite and does not need any regularization. 
Such a scheme is consistent with the usual regularization used in the solitonic calculations in the NJL model (see 
and the references therein) in which only the divergent part of the effective action is regularized, whereas the finite 
valence part is not affected by the regularization. Moreover, any regularization of the medium part would suppress 
the temperature effects ||, since the positive part of the spectrum would be affected by the cutoff as well. Thus 
regularizing only the divergent Dirac sea part, we obtain the effective action in the form 



S eff (fi,(3) = (3VN C [ [i# /a (ea) + (/* - O] + J E b [ 
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where the proper-time regularization function is given by 



R£(e) = 
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Since the Dirac sea does not contribute to the baryon number, expression (^) is finite and does not include any 
regularization. 
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III. FIXING THE MODEL PARAMETERS 



In the vacuum (/i = 0, T = 0) the stationary conditions (^|) lead to a translationary invariant solution 

a c = Mq and 7? c = , 



(16) 



which breaks the chiral symmetry. Due to the non-zero vacuum expectation value a c the quarks acquire a constituent 
mass Mq. The latter is a solution of the gap equation (stationary condition for the sigma field) which in the case of 
proper-time regularization has the form 



Nr. 



M = GMo—^? / 



2tt 2 



-2 



TO = mo - G < > 



(17) 



with < vp*? > being the quark condensate. Using (10) it is straightforward (see ref. |17| ) to evaluate the meson 
propagator: 
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The function Z p (Q 2 ) corresponds to a quark loop with two pseudoscalar-isovector insertions (375 7" Q ). Here we present 
only the proper-time regularized expression for it: 
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with r(0,x) being the incomplete gamma function. 

We fix the parameters of the model, namely the current mass too, the cutoff A and the coupling constant G in the 
vacuum sector reproducing the physical pion mass = 140 MeV and the pion decay constant f„ = 93 MeV. In fact, 
it leads to the Goldberger-Treiman (GT) relation on the quark level 



M = SW/tt , 

and one also recovers the Gell-Mann - Oakes - Renner (GMOR) relation: 

-too < > +0(toq) . 
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As usual, the last model parameter, the coupling constant G, is eliminated in favor of the constituent mass Mq usin 
the gap equation (|l7|). In fact, a value around Mo = 420 MeV is needed to describe properly the nucleon properties 
and we will use this value in our computations. The corresponding value of the proper-time cutoff is A = 640 MeV 



IV. PHASE TRANSITION AND MESON PROPERTIES AT FINITE TEMPERATURE AND DENSITY 



In this section we investigate the translationary invariant medium solution at finite temperature T and density p. 
For the model parameters mo, A and Mq we use the values fixed in the vacuum. 



A. Quark medium 

In this subsection we consider the medium as a Fermi sea of quarks. Solving (|^) together with the number of 
particles N in the volume V (0) kept fixed, we get the the medium solution in the form 

(7 c = M and tt c = . (22) 

As in the vacuum case, only the "classical" value of the scalar field can be non-zero and we call it the constituent 
quark mass M in medium. This mass M is a solution of the equation of motion the scalar field 
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FIG. 1. Constituent quark mass M as a function of temperature for different densities. At vanishing both temperature and 
density it is fixed to Mo = 420 MeV. The results with plane-wave basis (dotted lines) are compared to those of the finite-discrete 
basis (dashed lines). 
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with the proper-time regulator R^ 2 (e n ) in the form (|15|). The thermodynamical potential £1 is given by (| 
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Here, we subtract the vacuum value of the thermodynamical potential fl(0, 0). The energies are the eigenvalues of 
the hamiltonian h for the vacuum meson configuration (^|). Using (ffl) the baryon density pb can be written as 



1 3!) P 1 v ^ 
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and the free energy (||) is given by 

F{ii,0) = n(^l3)+N ci x PB . 



(25) 
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In the calculations we use a plane wave basis as well as a quasi-discrete basis and numerical method of Ripka and 
Kahana [^0| for solving the eigenvalue problem (||) by putting the system in a spherical box of a large radius D. The 
basis is made discrete by imposing a boundary condition at D. It is also made finite by introducing a numerical cut-off 
^max for the momenta of the basis states. Obviously both parameters are technical and the results do not depend on 
their particular values. The typical values which we use, are D « 18/Mo and k max w 8Mq. The volume V is taken to 
be a bit smaller than the box in order to avoid finite-size-box effects and the baryon number is given by 



N B = PB V. 



(27) 



The results for the constituent mass M as a function of temperature for different densities are presented in fig. [l] 
for both the plane wave and the quasi-discrete basis. As can be seen the quasi-discrete basis provides a quite good 
description of the continuum. 
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FIG. 2. Constituent quark mass M as a function of density for different temperatures. 

At finite density the constituent mass (or the chiral oder parameter < >) is a non-monotonic function of 
temperature. It means that in a hot matter M is less affected than in the case of cold matter. This also can be seen 
on fig. |^ where M is presented as function of density for different temperatures. At some critical values of temperature 
and/or density the constituent mass M is reduced to the current mass mo (the quark condensate < 'F'l' > vanishes) 
which is an indication for a a transition to the Wigner phase where the chiral symmetry is not spontaneously broken. 
The corresponding equation of state (pressure versus density) can be found in ref . |^] . The critical p — T chiral phase 
transition diagram is shown in fig. [|. At lower temperature values T < 90 MeV we have a clear first order transition 
which change to a second order one at higher temperatures. The critical value of the temperature is close to 200 MeV 
and for the density is about 2.5 p nm (nuclear matter density p nm — 0.16 fm~ 3 ). 

The meson masses in medium arc defined as the lowest zero solution of the inverse meson propagators 

at q = 0. The latter is given by the second variation of the effective action with respect to the meson fields 
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In the derivation of ([29|) we use the plane wave basis and eq.(£22|). V means principle value. From (|29j) it is clear 
that in the Goldstone phase the scalar and pseudoscalar mesons differ in mass. In the chiral limit the pions remain 
massless Goldstone bosons also in the medium, whereas the sigma mass follows the constituent quark mass. Both 
GMOR and GT relations are also valid in medium. Since in the Wigner phase the constituent mass M vanishes up 
to the current mass mo, the mesons become degenerated in mass and appear as parity-doubled mesons. In fact, in 
the Wigner phase due to the Fermi sea contribution, even in the chiral limit mo = the mesons have a non-zero 
mass, which increases with temperature (density). It is illustrated in fig. |I] where our results for the meson masses in 
medium are shown as a function of temperature at three different densities. Similar results are obtained by Zhuang 
et al. ]T^| in the case of finite temperature and vanishing chemical potential. 
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FIG. 3. p — T critical phase diagram. The lines separates the Goldstone phase with the chiral symmetry breaking from the 
Wigner phase where the chiral symmetry is restored. The solid line shows the critical values for which there is a first order 
transition, whereas the dashed line corresponds to the second order one. The arrow shows the temperature at which the order 
is changed. 



Since the NJL model lacks confinement there is a coupling to the qq continuum. In the Goldstone phase this 
coupling is unphysical but it does not affect much the structure of the scalar and pseudoscalar mesons - for the pion 
we have a <5-peak far from a very weak qq continuum and for the sigma meson there is a very sharp peak at the qq 
threshold (4M 2 ). The reason is that the pions play the role of Goldstone bosons and the sigma mesons are their chiral 
partner. In the Wigner phase, however, after the dcconfincment transition (suggested by the lattice results) the qq 
channel is physically open. The meson states are not more stable. 
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FIG. 4. Meson masses in the quark medium as a function of temperature for different densities. The shadowed area shows 
the width of the meson resonances. 
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They appear as resonances in the qq continuum and the poles of the meson propagators are complex. In this 
case one can identify the meson resonances with the peaks in the spectral density (imaginary part of the meson 
propagator). The spectral density is shown on fig. || for some different temperatures above the critical one. As can 
be seen, there is indeed a clear resonance structure. Close to the critical T cr as 200 MeV the peak is rather narrow. 
With the increasing temperature the resonance structures become weaker and at around 300 MeV is almost washed 
out. The position of the peak moves to higher energies and its width increases as well which implies that the mesons 
in the Wigner phase are rather soft modes than elementary excitations || . It is also interesting to notice that in this 
phase the meson masses show almost no density dependence. To conclude, according to the above model picture, at 
higher temperatures one might expect a smooth dehadronization transition. In this range, not the hadron but the 
quark and gluon degrees of freedom are expected to be relevant. Apparently, the missing explicit gluonic degrees of 
freedom in the NJL model limits its applicability at high temperatures. 




q 2 [GeV 2 ] 

FIG. 5. Imaginary part of the pion propagator (spectral density) for nuclear matter density and different temperatures. 

B. Medium of nucleons 

In the previous subsection we have considered the quark medium. In this case at finite temperature the single 
quarks from the Dirac sea are allowed be excited and occupy the levels in positive part of the spectrum (Fermi sea) 
leaving holes (antiquarks) in the Dirac sea. Apparently, because of confinement, at T and p below the critical values 
for the chiral and deconfinement transitions it is forbidden. As an alternative to this picture in this subsection we use 
a hybrid model in which instead of quarks we consider a Fermi sea of nucleons. The Dirac sea consists of quarks and 
it determines the vacuum sector. In this model the mesons are still qq excitations but they are also directly coupled 
to the nucleons of the Fermi sea. Such a model have been proposed and used by Jaminon et al. [fi~7| to study the 
the chiral symmetry restoration in a cold nucleon medium. At finite temperature only nucleons (colorless N c quark 
clusters) can be removed from the Dirac sea and occupy the levels in the Fermi sea. In this model picture, removing 
nucleons from the Dirac to the Fermi sea, one creates antinucleons. Apparently, it is physically reasonable to apply 
this hybrid model only for temperature and density below their critical values. As we see later, almost simultaneously 
with the chiral phase transition there is a delocalization of the soliton in the model, which means that after the phase 
transition it is consistent to consider a quark medium rather than the nucleon one. 

Since in the present scheme the contribution of the Dirac sea is separated from the one of the Fermi sea, it is 
straightforward to write down the medium contribution to the thermodynamical potential (effective action) in terms 
of nucleons: 
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FIG. 6. Constituent quark mass M as a function of temperature for different densities in the quark (dashed lines) and 
nucleon medium (solid lines). 
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The energies e„ are the solutions of the Dirac equation 
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The meson fields are coupled to the nucleons with a coupling constant <7at which relates the nucleon mass to the 
nonzero expectation value of the scalar meson field (constituent quarks mass Mq) in vacuum 



M N = g N M . 



(32) 



As in the quark spectrum, there is a gap of 2Mjv in the nucleon spectrum which separates the negative part of the 
spectrum from the positive one. In the nucleon medium, we also assume that g^ remains unchanged which means 
that a relation similar to ([32]) is valid also in medium. Since we are able to calculate the mass of the nucleon as a 
soliton, later we will check this approximation. The chemical potential fi is fixed by the baryon density 
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(or equivalently by the baryon number (|27D). It is interesting to make a parallel with the approach of Walecka et al. j25[| . 
Both models contain nucleons interacting locally with meson fields. However, in contrast to the Walecka model, where 
the meson fields are fundamental and the meson properties are introduced as phcnomcnological parameters, in the 
hybrid model the meson fields appears basically as qq pairs. Hence, we are able to calculate the meson masses and 
the physical coupling constants for a given constituent mass in vacuum as well as in medium. In the latter case, there 
is a change of the chiral condensate which leads to a modification of the meson properties as well. Moreover, within 
the present model the nucleon can be described as a soliton, which allows for a check of the basic assumption of the 
model, namely of the universality of the local coupling of the nucleons to the meson fields (^) . 

As in the case of a quark medium, only the scalar field can have a non-zero expectation value in nucleon medium 
(constituent quark mass in medium) which is a solution of the corresponding equation of motion ( p3|) with a finite- 
temperature part written now in terms of nucleons 
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FIG. 7. Constituent quark mass M as a function of density at T — in the quark (dashed line) and nucleon medium (solid 
line) . 
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together with the constraint pb = const. 

The results for M are presented in figj^ at finite density as a function of T and on fig.0 as a function of density 
at T = 0. At low T values the curves are close to those of the quark matter which means that the use of a quark 
matter instead of nucleon one is a reasonable approximation which is not true for higher temperatures. Already at 
intermediate temperatures T > 120 MeV they start to deviate significantly. Whereas the quark matter is quite soft 
against thermal fluctuations, the nucleon matter is much stiffer and the corresponding chiral phase transition is rather 
sharp. Both the chiral condensate and the constituent quark mass show a discontinuity at the critical temperature. 
In fact, the system jumps between two minima. The later is suggests that we have to deal with a first order phase 
transition even in the case of vanishing density in contrast to the case of quark matter where at temperatures T > 80 
MeV it is a second order one. The critical temperature is slightly larger than in the quark case and as we will see later 
it almost coincides with the critical temperature for the derealization of the soliton which makes this model picture 
consistent. The results at T = and finite p suggest a different model scenario. The critical density (w two times 
the nuclear matter density) for the chiral transition, which in fact, as we will see later, is also the critical density for 
the derealization of the soliton in cold nucleon matter, is smaller than the one for the quarks matter. This model 
scenario suggests a derealization transition from cold nucleon matter to quark one before the chiral phase transition 
to take place. 

According to the present model picture, at some critical temperature one expects a chiral phase transition together 
with a derealization transition from nucleon matter to quark one. Apparently, it would lead to drastic structural 
changes in the system. It can be clearly seen from the corresponding EOS (pressure versus baryon density) for 
different temperatures plotted in figj^. On this figure we combine the results from the nucleon matter (below the 
critical temperatures) in the hybrid model with those of the quark matter after the transition. Since we do not include 
vector mesons in the hybrid model, we are not able to reproduce the nuclear matter saturation at zero temperature 
and finite density which in the Walecka approach is due to the interplay between the cr-meson attraction and the 
u-meson repulsion. All curves in fig.|| show rapid change discontinuity at the derealization transitions from nucleon 
to quark matter. Similar behavior can be found on the fig. ^| and fig. [l^ for the pressure and energy density as a 
function of T for different densities. 

Apart from the medium part written in terms of nucleons ([34j) the inverse meson propagators in the nucleon medium 
have the same structure ( p9[ ) as in the quark medium. In chiral limit the pions are massless Goldstone bosons whereas 
the sigma mass is given by 2M. Like the constituent mass M the mass of the sigma meson in the nucleon medium 
has a discontinuity at the critical temperature. 
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FIG. 8. EOS (pressure versus density) for different temperatures. The pressure of the vacuum is subtracted. The disconti- 
nuities correspond to a phase transition from a nucleon to a quark matter. 
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FIG. 9. Pressure as a function of temperature for different densities. The discontinuities show the phase transition from the 
nucleon to quark matter. 

Our considerations are done in the mean-field approximation and suffer from the fact that meson quantum (loop) 
effects are not taken into account. However, as it is estimated in [Tl| , p^ ] the meson quantum effects may play a 
dominant role at low temperatures and vanishing density, whereas close to the critical temperature their contribution 
to the bulk thermodynamical variables is of order of 10 % or even less. Hence, following the result of |]ll|,^2), one 
expects that the meson fluctuations would not change the model scenario for the phase transition presented above. 



11 



> 

O 



I 

o 



10 



4 - 



2 - 



Mr, = 420 McV 



p = Q 

p = 1 Pnn 
P = 2 pun 



:7 



l l I l l l l I l l l l_ 



' I l l l l I l l l l_ 



■30 



100 150 200 

Temperature [MeV] 



250 



300 



FIG. 10. Energy density as a function of temperature for different densities. The discontinuities show the phase transition 
from the nucleon to quark matter. 



V. NUCLEON AS A SOLITON IN MEDIUM 



As a next step we solve the B = 1 soliton sector in a hot medium. We look for a localized bound solution (soliton) 
of N c valence quarks interacting with the Fermi and Dirac sea both polarized due to the interaction. 



A. In quark medium 



First we consider the case of Fermi and Dirac sea both filled with quarks of a constituent mass M. The thermody- 
namical potential (effective action) includes an explicit valence quark contribution: 



(3) = N c 6(e val )e val + N c £ [i^ /2 (e„) - R$ /2 (e° n ) + (/i 

e„<0 



1 

V 



^(M 2 -M 2 )-^(a-M ) 



since the valence level is always occupied by N c quarks. 
The meson fields are assumed to be in a hedgehog form 



T ^ 

e n ^val 



l 



-/3(e„-p) 



(35) 



a(r) = <j{r) 



and 



7?(r) = r7r(r) 



(36) 



and arc restricted on the chiral circle 



(7 + 7T 



M 



(37) 



The latter is an ad hoc nonlinear constraint which reflects the fact that the pions as almost massless Goldstone bosons 
play a privileged role of dynamical mesons in this low-energy region pjj . 
From (||) one gets the equations of motion for the meson fields 
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FIG. 11. Soliton energy as a function of the temperature at three different quark matter densities. 
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(39) 



We use a numerical self-consistent iterative procedure based on a method proposed by Ripka and Kahana [g0| . The 
procedure consists in solving in an iterative way the Dirac equation (||) together with the equations of motion of the 
meson fields - eqs.(|3"8|),(|39|), and the constraint Nb = const (]27]) which is actually a condition fixing the chemical 
potential /i. 

The techniques for the numerical procedure are well known |^,^3| for the case of vanishing chemical potential 
where there is no contribution coming from the positive continuum. They can be easily adopted to the case of finite 
chemical potential and temperature where an additional source (Fermi sea quarks) for the meson fields appears in 
eqs.dH),©. 

In the case of fixed T and p the proper way to describe the equilibrium state of a thermodynamical system is to use 
the free energy (^). Hence, the energy of the B — 1 soliton (effective soliton mass) is given by the change of the free 
energy when the N c valence quarks are added to the medium. Apparently, in this the baryon number of the system 
increases by one, Nb + 1- Since due to the presence of the soliton both the Fermi sea and the Dirac sea are getting 
polarized, subtracting the free energy F(fiQ,{3) of the unperturbed Fermi and Dirac sea (translationally invariant 
medium solution), the soliton energy is given by the sum of the energy of the valence quarks and the contributions 
due to the polarization of both continua: 



E so i = N c ri va ie va i + N c 
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FIG. 12. Soliton energy as a function of density for vanishing as well as finite temperature values. 

The chemical potential u fixed by the constrain Ng — const ( p7| ) corresponds to the translationally invariant 
medium solution. The soliton baryon density distribution is split in valence, sea and medium parts: 



Ps 



oi(r) = 9(e val )<S>l al (r)<P val (r) + < (f)$„(r)sgn(-e n ) 



2-^1 l _|_ e (e„-/i)/3 1 . 

t n >M e„<0 



(41) 



Because of the constraint Nb = const, only the first term in r.h.s. of ( [4l| ) contributes to the baryon number B = 1. 
The soliton m.s. radius defined as 



< r 2 > so i= J d 3 rr 2 p B =i{r) 



(42) 



can be used to measure the spatial extension of the soliton. 

Our results for the soliton energy as a function of the temperature for vanishing as well as for finite density values 
are presented in fig. [ll]. For completeness on the fig. [l^ we also present the soliton energy as a function of density at 
vanishing [ ]l4[ as well as at finite temperature. 

We do not find a localized solution (soliton) at temperatures larger than a critical value Tf r ss 180 MeV. It means 
that at large enough values the temperature effects simply disorder the system and destroys the soliton. Within the 
present model picture this may be interpreted as an indication for a derealization of the nucleon in hot medium but, 
however, one should keep in mind that the model lacks confinement. It should be mentioned that a similar effect has 
been found in the Skyrme model |24j . 

As can be seen from fig. |ll|the temperature effects for the soliton are much weaker than the finite density effects. In 
the case of finite T and zero p the soliton energy shows almost no reduction for temperatures not close to the critical 
one, whereas in the opposite case of finite p and zero T the soliton energy is linearly decreasing with the density. In 
the hot medium (both T and p finite) the temperature clearly suppresses the finite density effects and stabilizes the 
soliton. The latter results in a non-monotonic temperature dependence of the soliton energy at fixed finite density. 
At densities larger than the nuclear matter one p nm , which are relevant for the heavy-ion experiments, the reduction 
of the soliton energy at finite temperature is much smaller than one at T — (see fig. O). 
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FIG. 13. Soliton energy and contributions coming from the valence quarks and from the polarized continua as a function of 
temperature at density p = p nm - 

In particular, at p nm and T > 100 MeV the reduction is about 15% whereas at T = is two times larger. Further, 
at densities larger than two times p nm the soliton exists only at intermediate temperatures 100 MeV < T < 180 MeV. 
It means that according to the present model calculations, the soliton is more stable in the hot matter than in the 
cold one. In fact, it can be easily understood. Because of the gap in the quark spectrum, 2M, the Dirac sea is much 
less affected by the temperature than the Fermi sea (positive-energy part of the spectrum). Disordering mainly the 
Fermi sea the attractive interaction between the medium and the valence quarks and Dirac sea, which destabilize the 
soliton, is diminished. 
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At temperatures close to the critical one, the thermal fluctuations become comparable with the chiral order param- 
eter, the chiral condensate qq >, they completely disorder the system and in particular, destroy the soliton. 

It is interesting to look at the separated valence and sea contributions to the soliton energy presented in fig. [l3| 
and [H] for two different medium densities, namely one and two times nuclear matter density, respectively. Whereas 
at low temperature almost a half of the energy comes from the polarized Dirac and Fermi sea, at increasing T the 
valence contribution becomes dominant. Also, in contrast to the valence quark contribution, which shows a strong 
dependence on T, the polarized sea contribution stays almost constant with T not close to the critical one but it 
is strongly affected by the medium density. In particular, going from one to two times p n m, the sea contribution is 
reduced by a factor of two, which destabilizes the soliton. Further, at T close to the critical values it vanishes. 

In order to illustrate the change of the soliton structure in hot medium we also plot the soliton square radius (^2|) as 
a function of density for different temperatures in figs. [l5|. All curves in fig. [l^ show a clear trend to grow rapidly at 
temperatures close to the critical values which is an indication for derealization of the soliton. At finite both density 
and temperature the radius is smaller than in the case of cold matter which is a sign for a stabilization of the soliton 
in hot medium compared to the case of cold one. At densities larger than two times p nm the soliton exists only at 
intermediate temperatures. 
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FIG. 15. Soliton m.s. radius as a function of density for vanishing as well as for finite temperature values. 

B. In medium of nucleons 

In this subsection we consider the B = 1 soliton in the hybrid model with N c valence quarks interacting with both 
the quark Dirac sea and the nucleon Fermi sea via the meson fields. As in the quark medium the B — 1 soliton is 
obtained solving the Dirac equations (|^),(^l|) together with meson equations of motion (^8|) , (|39|) which contain the 
finite-temperature nucleon part (|34|), and the constraint pb = const ([33]) to fix the chemical potential p N in the 
self-consistent iterative procedure. The energy of the B = 1 soliton is given again by the change of the free energy 
when the N c valence quarks are added to the system. 

On fig.|l6| the temperature dependence of the calculated B = 1 soliton energy for the nuclear matter density is 
compared with those of the quark matter. The two curves have similar trends at intermediate values of T and start 
to deviate at T close to critical one. The soliton in the nucleon matter is more bound and less affected by the 
temperature. At some critical temperature of about 200 MeV the soliton disappears - derealization of the soliton. In 
fact, at the same temperature in the hybrid model there is also the chiral phase transition from Goldstone to Wigner 
phase. The m.s. soliton radius is presented in fig.[l7]. Close to the critical temperature it starts to grow but it is much 
less pronounced than in the case of quark matter. 
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FIG. 16. Soliton energy in a nucleon matter compared to those in a quark matter for p = p nm . 

Back to our discussion concerning the Walecka model, in our case this model picture is valid to some critical density 
of about two times nuclear matter density. At higher densities we found a derealization transition to a quark matter. 
On hg.ttq we also plot the coupling constant 



g N = E soi /M 



(43) 



as a function of T at one nuclear matter density. As can be seen it stays almost constant which means that the 
relation (32) is a good approximation also in hot medium. 
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In contrast to the quark matter the energy of the soliton in the nucleon matter is slightly above the threshold of 
3Af . However, in the hybrid model this channel is physically closed by construction - before the phase transition, in 
the Goldstone phase, free quarks are not allowed to occupy the positive part of the quark spectrum. 
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FIG. 18. Ratio E ao \/M in a nucleonic medium for p — p nm as a function of temperature. 



VI. SUMMARY 



We study the bulk thermodynamical properties, some meson properties and the nucleon as a B = 1 soliton in hot 
medium using the NJL model in mean-field approximation. We consider the case of a quark medium in the NJL model 
as well as of a nucleon medium in a hybrid NJL model in which the quark Dirac sea is combined with a nucleon Fermi 
sea. In both cases we find that the chiral order parameter, the chiral condensate, gets modified and at some critical 
values of temperature and/or density vanishes which indicates a chiral phase transition from Goldstone to Wigner 
phase. At finite density the chiral order parameter and the constituent quark mass have a non-monotonic temperature 
dependence and at finite temperatures not close to the critical one they are less affected than in the cold matter. 
Below the chiral phase transition the pions are Goldstone bosons and their mass remains almost unchanged whereas 
the sigma mass follows the constituent quark mass. After the phase transition the mesons are parity-doubled and 
appear as broad resonances in the qq continuum. With increasing temperature the position of the resonance moves 
to higher energies and and its width increases very fast almost independently of the density. At some temperature of 
about 300 MeV the meson structure is very weak and is almost washed out. Whereas at low density and temperature 
values the quark matter and the nucleon matter provide similar results, at larger temperatures they differ significantly. 
In particular, the quark matter is much softer against thermal fluctuations and the chiral phase transition is rather 
smooth. The nucleon matter is much less affected by the temperatures even for values close to the critical one. The 
chiral phase transition is rather sharp and all thermodynamical variables show a large discontinuities which is an 
indication for a first order phase transition. We study also the structure of the baryon number one soliton of the NJL 
model immersed in a hot quark as well as in a nucleon medium. The polarization of both the Fermi and the Dirac 
sea is taken into account in a consistent way. We find that at finite density the temperature stabilizes the soliton 
and after some critical density it exists only at intermediate temperature values. In general, at finite temperature 
the soliton is more bound and less swelled than in the case of a cold matter. At some critical temperature we find 
no more a localized solution. In the case of nucleon matter this derealization means a transition from the nucleon 
to the quark matter. The critical temperature coincides with those for the chiral phase transition. According to this 
model scenario one should expects at some critical temperature a common sharp phase transition from the nucleon 
to the quark matter with a restoration of chiral symmetry. All present results are obtained in mean-field (lar ge N c ) 
approximation which means that the meson quantum (loop) effects are not included. However, it is expected J11 12 
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that these effects would play a minor role near the phase transition and hence it would not change principally our 
results. 
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